A re-activation model for 169Yb intensity modulated brachytherapy sources accounting for spatiotemporal isotopic composition.

BACKGROUND: Intensity modulated brachytherapy based on partially shielded intracavitary and interstitial applicators is possible with a cost-effective 169Yb production method. 169Yb is a traditionally expensive isotope suitable for this purpose, with an average γ-ray energy of 93 keV. Re-activating a single 169Yb source multiple times in a nuclear reactor between clinical uses was shown to theoretically reduce cost by approximately 75% relative to conventional single-activation sources. With re-activation, substantial spatiotemporal variation in isotopic source composition is expected between activations via 168Yb burnup and 169Yb decay, resulting in time dependent neutron transmission, precursor usage, and reactor time needed per re-activation. PURPOSE: To introduce a generalized model of radioactive source production that accounts for spatiotemporal variation in isotopic source composition to improve the efficiency estimate of the 169Yb production process, with and without re-activation. METHODS AND MATERIALS: A time-dependent thermal neutron transport, isotope transmutation, and decay model was developed. Thermal neutron flux within partitioned sub-volumes of a cylindrical active source was calculated by raytracing through the spatiotemporal dependent isotopic composition throughout the source, accounting for thermal neutron attenuation along each ray. The model was benchmarked, generalized, and applied to a variety of active source dimensions with radii ranging from 0.4 to 1.0 mm, lengths from 2.5 to 10.5 mm, and volumes from 0.31 to 7.85 mm3, at thermal neutron fluxes from 1 × 1014 to 1 × 1015 n cm-2 s-1. The 168Yb-Yb2O3 density was 8.5 g cm-3 with 82% 168Yb-enrichment. As an example, a reference re-activatable 169Yb active source (RRS) constructed of 82%-enriched 168Yb-Yb2O3 precursor was modeled, with 0.6 mm diameter, 10.5 mm length, 3 mm3 volume, 8.5 g cm-3 density, and a thermal neutron activation flux of 4 × 1014 neutrons cm-2 s-1. RESULTS: The average clinical 169Yb activity for a 0.99 versus 0.31 mm3 source dropped from 20.1 to 7.5 Ci for a 4 × 1014 n cm-2 s-1 activation flux and from 20.9 to 8.7 Ci for a 1 × 1015 n cm-2 s-1 activation flux. For thermal neutron fluxes ≥2 × 1014 n cm-2 s-1, total precursor and reactor time per clinic-year were maximized at a source volume of 0.99 mm3 and reached a near minimum at 3 mm3. When the spatiotemporal isotopic composition effect was accounted for, average thermal neutron transmission increased over RRS lifetime from 23.6% to 55.9%. A 28% reduction (42.5 days to 30.6 days) in the reactor time needed per clinic-year for the RRS is predicted relative to a model that does not account for spatiotemporal isotopic composition effects. CONCLUSIONS: Accounting for spatiotemporal isotopic composition effects within the RRS results in a 28% reduction in the reactor time per clinic-year relative to the case in which such changes are not accounted for. Smaller volume sources had a disadvantage in that average clinical 169Yb activity decreased substantially below 20 Ci for source volumes under 1 mm3. Increasing source volume above 3 mm3 adds little value in precursor and reactor time savings and has a geometric disadvantage.

Integration and dosimetric validation of a dynamic collimation system for pencil beam scanning proton therapy.

Objective.To integrate a Dynamic Collimation System (DCS) into a pencil beam scanning (PBS) proton therapy system and validate its dosimetric impact.Approach.Uncollimated and collimated treatment fields were developed for clinically relevant targets using an in-house treatment plan optimizer and an experimentally validated Monte Carlo model of the DCS and IBA dedicated nozzle (DN) system. The dose reduction induced by the DCS was quantified by calculating the mean dose in 10- and 30-mm two-dimensional rinds surrounding the target. A select number of plans were then used to experimentally validate the mechanical integration of the DCS and beam scanning controller system through measurements with the MatriXX-PT ionization chamber array and EBT3 film. Absolute doses were verified at the central axis at various depths using the IBA MatriXX-PT and PPC05 ionization chamber.Main results.Simulations demonstrated a maximum mean dose reduction of 12% for the 10 mm rind region and 45% for the 30 mm rind region when utilizing the DCS. Excellent agreement was observed between Monte Carlo simulations, EBT3 film, and MatriXX-PT measurements, with gamma pass rates exceeding 94.9% for all tested plans at the 3%/2 mm criterion. Absolute central axis doses showed an average verification difference of 1.4% between Monte Carlo and MatriXX-PT/PPC05 measurements.Significance.We have successfully dosimetrically validated the delivery of dynamically collimated proton therapy for clinically relevant delivery patterns and dose distributions with the DCS. Monte Carlo simulations were employed to assess dose reductions and treatment planning considerations associated with the DCS.

The population percentile allowance method for determining systematic spatial error tolerances for temporary intensity modulated brachytherapy.

BACKGROUND: Multiple approaches are under development for delivering temporary intensity modulated brachytherapy (IMBT) using partially shielded applicators wherein the delivered dose distributions are sensitive to spatial uncertainties in both the applicator position and shield orientation, rather than only applicator position as with conventional high-dose-rate brachytherapy (HDR-BT). Sensitivity analyses to spatial uncertainties have been reported as components of publications on these emerging technologies, however, a generalized framework for the rigorous determination of the spatial uncertainty tolerances of dose-volume parameters is needed. PURPOSE: To derive and present the population percentile allowance (PPA) method, a generalized mathematical and statistical framework to evaluate the tolerance of temporary IMBT approaches to spatial uncertainties in applicator position and shield orientation. METHODS: A mathematical formalism describing geometric applicator position and shield orientation shifts was derived that supports straight and curved applicators and applies to serial and helical rotating shield brachytherapy (RSBT) and direction modulated brachytherapy (DMBT). The PPA method entails defining the percentage of a patient population receiving a given therapy that is, allowed to receive dose-volume errors in the target volume and specified organs at risk of a defined percentage or less, then determining what combinations of applicator position and shield orientation systematic errors would be expected to produce that outcome in the population. The PPA method was applied to the use case of multi-shield helical 169 Yb-based RSBT for cervical cancer, with 45° and 180° shield emission angles. A total of 37 cervical cancer patients were considered in the population, with average (± 1 standard deviation) HR-CTV volumes of 79 cm3  ± 37 cm3 and optimized baseline treatment plans (no spatial uncertainties applied) created for each patient to meet dose-volume requirements of 85 GyEQD2 (equivalent uniform dose in 2 Gy fraction), with D2cc tolerance doses of 90 GyEQD2 , 75 GyEQD2 , and 75 GyEQD2 for bladder, rectum, and sigmoid colon, respectively. RESULTS: For the PPA requirement that 90% of cervical cancer patients receiving multi-shield helical RSBT could have a maximum dose-volume uncertainty of 10% for high-risk clinical target volume (HR-CTV) D90 (minimum dose to hottest 90%) and bladder, rectum, and sigmoid colon D2cc (minimum dose to hottest 2 cm3 ), the tolerance systematic applicator position and shield orientation uncertainties were approximately ± 1.0 mm and ± 4.25°, respectively. For ± 1.5 mm and ± 5° systematic applicator position and shield orientation tolerances, 90% of the patients considered would have a maximum dose-volume uncertainty of 12.8% or less. CONCLUSION: The PPA method was formalized to determine the temporary IMBT spatial uncertainty tolerances that would be expected to result in an allowed percentage of a population of patients receiving relative dose-volume errors above a defined percentage. Multi-shield, helical 169 Yb-based RSBT for cervical cancer was evaluated and tolerances determined, which, if applied on each treatment fraction, would represent an extreme situation. The PPA method is applicable to a variety of temporary IMBT approaches and can be used to rigorously determine the design parameters for the delivery systems such as mechanical driver motor accuracy, shield angle backlash, applicator rotation, and applicator fixation stability.

Design, testing and characterization of a proton central axis alignment device for the dynamic collimation system.

Objective. Proton therapy conformity has improved over the years by evolving from passive scattering to spot scanning delivery technologies with smaller proton beam spot sizes. Ancillary collimation devices, such the Dynamic Collimation System (DCS), further improves high dose conformity by sharpening the lateral penumbra. However, as spot sizes are reduced, collimator positional errors play a significant impact on the dose distributions and hence accurate collimator to radiation field alignment is critical.Approach. The purpose of this work was to develop a system to align and verify coincidence between the center of the DCS and the proton beam central axis. The Central Axis Alignment Device (CAAD) is composed of a camera and scintillating screen-based beam characterization system. Within a light-tight box, a 12.3-megapixel camera monitors a P43/Gadox scintillating screen via a 45° first-surface mirror. When a collimator trimmer of the DCS is placed in the uncalibrated center of the field, the proton radiation beam continuously scans a 7×7 cm2square field across the scintillator and collimator trimmer while a 7 s exposure is acquired. From the relative positioning of the trimmer to the radiation field, the true center of the radiation field can be calculated.Main results.The CAAD can calculate the offset between the proton beam radiation central axis and the DCS central axis within 0.054 mm accuracy and 0.075 mm reproducibility.Significance.Using the CAAD, the DCS is now able to be aligned accurately to the proton radiation beam central axis and no longer relies on an x-ray source in the gantry head which is only validated to within 1.0 mm of the proton beam.

Dosimetric delivery validation of dynamically collimated pencil beam scanning proton therapy.

Objective. Pencil beam scanning (PBS) proton therapy target dose conformity can be improved with energy layer-specific collimation. One such collimator is the dynamic collimation system (DCS), which consists of four nickel trimmer blades that intercept the scanning beam as it approaches the lateral extent of the target. While the dosimetric benefits of the DCS have been demonstrated through computational treatment planning studies, there has yet to be experimental verification of these benefits for composite multi-energy layer fields. The objective of this work is to dosimetrically characterize and experimentally validate the delivery of dynamically collimated proton therapy with the DCS equipped to a clinical PBS system.Approach. Optimized single field, uniform dose treatment plans for 3 × 3 × 3 cm3target volumes were generated using Monte Carlo dose calculations with depths ranging from 5 to 15 cm, trimmer-to-surface distances ranging from 5 to 18.15 cm, with and without a 4 cm thick polyethylene range shifter. Treatment plans were then delivered to a water phantom using a prototype DCS and an IBA dedicated nozzle system and measured with a Zebra multilayer ionization chamber, a MatriXX PT ionization chamber array, and Gafchromic™ EBT3 film.Main results. For measurements made within the SOBPs, average 2D gamma pass rates exceeded 98.5% for the MatriXX PT and 96.5% for film at the 2%/2 mm criterion across all measured uncollimated and collimated plans, respectively. For verification of the penumbra width reduction with collimation, film agreed with Monte Carlo with differences within 0.3 mm on average compared to 0.9 mm for the MatriXX PT.Significance. We have experimentally verified the delivery of DCS-collimated fields using a clinical PBS system and commonly available dosimeters and have also identified potential weaknesses for dosimeters subject to steep dose gradients.

Mechanical Characterization and Validation of the Dynamic Collimation System Prototype for Proton Radiotherapy.

Radiation therapy is integral to cancer treatments for more than half of patients. Pencil beam scanning (PBS) proton therapy is the latest radiation therapy technology that uses a beam of protons that are magnetically steered and delivered to the tumor. One of the limiting factors of PBS accuracy is the beam cross-sectional size, similar to how a painter is only as accurate as the size of their brush allows. To address this, collimators can be used to shape the beam along the tumor edge to minimize the dose spread outside of the tumor. Under development is a dynamic collimation system (DCS) that uses two pairs of nickel trimmers that collimate the beam at the tumor periphery, limiting dose from spilling into healthy tissue. Herein, we establish the dosimetric and mechanical acceptance criteria for the DCS based on a functioning prototype and Monte Carlo methods, characterize the mechanical accuracy of the prototype, and validate that the acceptance criteria are met. From Monte Carlo simulations, we found that the trimmers must be positioned within ±0.5 mm and ±1.0 deg for the dose distributions to pass our gamma analysis. We characterized the trimmer positioners at jerk values up to 400 m/s3 and validated their accuracy to 50 µm. We measured and validated the rotational trimmer accuracy to ±0.5 deg with a FARO® ScanArm. Lastly, we calculated time penalties associated with the DCS and found that the additional time required to treat one field using the DCS varied from 25-52 s.

The dosimetric enhancement of GRID profiles using an external collimator in pencil beam scanning proton therapy.

PURPOSE: The radiobiological benefits afforded by spatially fractionated (GRID) radiation therapy pairs well with the dosimetric advantages of proton therapy. Inspired by the emergence of energy-layer specific collimators in pencil beam scanning (PBS), this work investigates how the spot spacing and collimation can be optimized to maximize the therapeutic gains of a GRID treatment while demonstrating the integration of a dynamic collimation system (DCS) within a commercial beamline to deliver GRID treatments and experimentally benchmark Monte Carlo calculation methods. METHODS: GRID profiles were experimentally benchmarked using a clinical DCS prototype that was mounted to the nozzle of the IBA-dedicated nozzle system. Integral depth dose (IDD) curves and lateral profiles were measured for uncollimated and GRID-collimated beamlets. A library of collimated GRID dose distributions were simulated by placing beamlets within a specified uniform grid and weighting the beamlets to achieve a volume-averaged tumor cell survival equivalent to an open field delivery. The healthy tissue sparing afforded by the GRID distribution was then estimated across a range of spot spacings and collimation widths, which were later optimized based on the radiosensitivity of the tumor cell line and the nominal spot size of the PBS system. This was accomplished by using validated models of the IBA universal and dedicated nozzles. RESULTS: Excellent agreement was observed between the measured and simulated profiles. The IDDs matched above 98.7% when analyzed using a 1%/1-mm gamma criterion with some minor deviation observed near the Bragg peak for higher beamlet energies. Lateral profile distributions predicted using Monte Carlo methods agreed well with the measured profiles; a gamma passing rate of 95% or higher was observed for all in-depth profiles examined using a 3%/2-mm criteria. Additional collimation was shown to improve PBS GRID treatments by sharpening the lateral penumbra of the beamlets but creates a trade-off between enhancing the valley-to-peak ratio of the GRID delivery and the dose-volume effect. The optimal collimation width and spot spacing changed as a function of the tumor cell radiosensitivity, dose, and spot size. In general, a spot spacing below 2.0 cm with a collimation less than 1.0 cm provided a superior dose distribution among the specific cases studied. CONCLUSIONS: The ability to customize a GRID dose distribution using different collimation sizes and spot spacings is a useful advantage, especially to maximize the overall therapeutic benefit. In this regard, the capabilities of the DCS, and perhaps alternative dynamic collimators, can be used to enhance GRID treatments. Physical dose models calculated using Monte Carlo methods were experimentally benchmarked in water and were found to accurately predict the respective dose distributions of uncollimated and DCS-collimated GRID profiles.

Development and validation of the Dynamic Collimation Monte Carlo simulation package for pencil beam scanning proton therapy.

PURPOSE: The aim of this work was to develop and experimentally validate a Dynamic Collimation Monte Carlo (DCMC) simulation package specifically designed for the simulation of collimators in pencil beam scanning proton therapy (PBS-PT). The DCMC package was developed using the TOPAS Monte Carlo platform and consists of a generalized PBS source model and collimator component extensions. METHODS: A divergent point-source model of the IBA dedicated nozzle (DN) at the Miami Cancer Institute (MCI) was created and validated against on-axis commissioning measurements taken at MCI. The beamline optics were mathematically incorporated into the source to model beamlet deflections in the X and Y directions at the respective magnet planes. Off-axis measurements taken at multiple planes in air were used to validate both the off-axis spot size and divergence of the source model. The DCS trimmers were modeled and incorporated as TOPAS geometry extensions that linearly translate and rotate about the bending magnets. To validate the collimator model, a series of integral depth dose (IDD) and lateral profile measurements were acquired at MCI and used to benchmark the DCMC performance for modeling both pristine and range shifted beamlets. The water equivalent thickness (WET) of the range shifter was determined by quantifying the shift in the depth of the 80% dose point distal to the Bragg peak between the range shifted and pristine uncollimated beams. RESULTS: A source model of the IBA DN system was successfully commissioned against on- and off-axis IDD and lateral profile measurements performed at MCI. The divergence of the source model was matched through an optimization of the source-to-axis distance and comparison against in-air spot profiles. The DCS model was then benchmarked against collimated IDD and in-air and in-phantom lateral profile measurements. Gamma analysis was used to evaluate the agreement between measured and simulated lateral profiles and IDDs with 1%/1 mm criteria and a 1% dose threshold. For the pristine collimated beams, the average 1%/1 mm gamma pass rates across all collimator configurations investigated were 99.8% for IDDs and 97.6% and 95.2% for in-air and in-phantom lateral profiles. All range shifted collimated IDDs passed at 100% while in-air and in-phantom lateral profiles had average pass rates of 99.1% and 99.8%, respectively. The measured and simulated WET of the polyethylene range shifter was determined to be 40.9 and 41.0 mm, respectively. CONCLUSIONS: We have developed a TOPAS-based Monte Carlo package for modeling collimators in PBS-PT. This package was then commissioned to model the IBA DN system and DCS located at MCI using both uncollimated and collimated measurements. Validation results demonstrate that the DCMC package can be used to accurately model other aspects of a DCS implementation via simulation.

Needle-free cervical cancer treatment using helical multishield intracavitary rotating shield brachytherapy with the 169 Yb Isotope.

PURPOSE: To assess the capability of an intracavitary 169 Yb-based helical multishield rotating shield brachytherapy (RSBT) delivery system to treat cervical cancer. The proposed RSBT delivery system contains a pair of 1.25 mm thick platinum partial shields with 45° and 180° emission angles, which travel in a helical pattern within the applicator. METHODS: A helically threaded tandem applicator with a 45° tandem curvature containing a helically threaded catheter was designed. A 0.6 mm diameter 169 Yb source with a length of 10.5 mm was simulated. A 37-patient treatment planning study, based on Monte Carlo dose calculations using MCNP5, was conducted with high-risk clinical target volumes (HR-CTVs) of 41.2-192.8 cm3 (average ± standard deviation of 79.9 ± 35.8 cm3 ). All patients were assumed to receive 25 fractions of 1.8 Gy of external beam radiation therapy (EBRT) before receiving 5 fractions of high-dose-rate brachytherapy (HDR-BT). For each patient, 192 Ir-based intracavitary (IC) HDR-BT, 192 Ir-based intracavitary/interstitial (IC/IS) HDR-BT using a hybrid applicator with eight IS needles, and 169 Yb-based RSBT plans were generated. RESULTS: For the IC, IC/IS, and RSBT treatment plans, 38%, 84%, and 86% of the plans, respectively, met the planning goal of an HR-CTV D90 (minimum dose to hottest 90%) of 85 GyEQD2 (α/ß = 10 Gy). Median (25th percentile, 75th percentile) treatment times for IC, IC/IS, and RSBT were 11.71 (6.62, 15.40) min, 68.00 (45.02, 80.02) min, and 25.30 (13.87, 35.39) min, respectively. 192 Ir activities ranging from 159.1-370 GBq (4.3-10 Ci) and 169 Yb activities ranging from 429.2-999 GBq (11.6-27 Ci) were used, which correspond to the same clinical ranges of dose rates at 1 cm off-source-axis in water. Extra needle insertion and planning time beyond that needed for intracavitary-only approaches was accounted for in the IC/IS treatment time calculations. CONCLUSION: 169 Yb-based RSBT for cervical cancer met the HR-CTV D90 goal of 85 Gy in a greater percentage of the patients considered than IC/IS (86% vs 84%, respectively) and can reduce overall treatment time relative to IC/IS. 169 Yb-based RSBT could be used to replace IC/IS in instances where IC/IS treatment is not available, especially in instances when HR-CTV volumes are ≥30 cm3 .

Development of a radiobiological evaluation tool to assess the expected clinical impacts of contouring accuracy between manual and semi-automated segmentation algorithms.

RADEval is a tool developed to assess the expected clinical impact of contouring accuracy when comparing manual contouring and semi-automated segmentation. The RADEval tool, designed to process large scale datasets, imported a total of 2,760 segmentation datasets, along with a Simultaneous Truth and Performance Level Estimation (STAPLE) to act as ground truth tumor segmentations. Virtual dose-maps were created within RADEval and two different tumor control probability (TCP) values using a Logistic and a Poisson TCP models were calculated in RADEval using each STAPLE and each dose-map. RADEval also virtually generated a ring of normal tissue. To evaluate clinical impact, two different uncomplicated TCP (UTCP) values were calculated in RADEval by using two TCP-NTCP correlation parameters (δ = 0 and 1). NTCP values showed that semi-automatic segmentation resulted in lower NTCP with an average 1.5 - 1.6 % regardless of STAPLE design. This was true even though each normal tissue was created from each STAPLE (p <; 0.00001). TCP and UTCP presented no statistically significant differences (p ≥ 0.1884). The intra-operator standard deviations (SDs) for TCP, NTCP and UTCP were significantly lower for the semi-automatic segmentation method regardless of STAPLE design (p <; 0.0331). Both intra-and inter-operator SDs of TCP, NTCP and UTCP were significantly lower for semi-automatic segmentation for the STAPLE 1 design (p <;0.0331). RADEval was able to efficiently process 4,920 datasets of two STAPLE designs and successfully assess the expected clinical impact of contouring accuracy.